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O Abstract 

o 

^vq We consider a quantum dot embedded in a three-dimensional nanowire with tunable aspect ratio a. A configuration interaction 
theory is developed to calculate the energy spectra of the finite ID quantum dot systems charged with two electrons in the presence 
of magnetic fields B along the wire axis. Fruitful singlet-triplet transition behaviors are revealed and explained in terms of the 
£^ competing exchange interaction, correlation interaction, and spin Zeeman energy. In the high aspect ratio regime, the singlet-triplet 
transitions are shown designable by tuning the parameters a and B. The transitions also manifest the highly correlated nature of 
long nanowire quantum dots. 



Key words: Nanowire quantum dot, Exchange, Correlation, Singlet-triplet transition 



i 



c3 



i 

C 

O 

o 



> 

00 

oo 
(N 

o 



X 



1. Introduction 

For years, few electron charged quantum dots have attracted 
extensive attention due to the controllable electronic and spin 
properties Jit]. However, only few attempts have so far 
been made for studies of finite ID nanowire quantum dots 
(NWQDs). More recently, it was shown that the NWQDs 
formed in the heterostructures in nanowires can be fabricated as 
single electron transistors and successively charged with con- 
trolled number of electrons [2]. The successful experimental 
works motivate us to explore the possible geometric effects of 
NWQDs characterized by their aspect ratios, a, on the elec- 
tronic and spin properties of two-electron charged NWQDs. 

In this work, we focus on the study of the singlet-triplet (ST) 
transitions in two-electron charged NWQDs yj], conducted by 
using a developed configuration interaction (CI) theory in com- 
bination with the exact diagonalization techniques based on a 
3D asymmetric parabolic model. It will be illustrated that the 
ST transitions in InAs-based NWQDs driven by an appropri- 
ate magnetic field are associated with the competing effects of 
large spin-Zeeman energies as well as the exchange and cor- 
relation energies. The correlation-dominated nature of a long 
NWQD (i.e. with high aspect ratio) will be identified by the 
spin phase diagram with respect to the applied magnetic fields 
and the tunable aspect ratio. 

2. Theoretical Model 

2.1. single-electron spectrum 

We begin with the single electron problem of a NWQD with 
axial magnetic field B = (0, 0, B), described by the Hamilto- 



nian J3l 



Ho = ^t(P + ehf + V(x,y,z) + H z , 
2m 



(1) 



where A = (B/2)(y, -x, 0) denotes the vector potential and m* 
stands for the effective mass of an electron with charge —e. 
The spin-Zeeman Hamiltonian Hz = g*/UBBs z is m terms of 
the z-component of electron spin s z and the effective Lande g- 
factor of electron g* and the Bohr magneton /ig. In addition, 
the confining potential V(x,y,z) - m* |wq (x 2 + y 2 ) + w 2 z 2 ] /2 
is assumed of the parabolic form with ojq and a> z parametriz- 
ing, respectively, the transverse and the longitudinal confining 
strength. 

In this work, we assume a constant g* and take the g* = -8.0 
for InAs The single electron Hamiltonian (Q3 leads to the 

extended Fock-Darwin single-particle spectrum 
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where n,m,q = 0, 1 , 2 ■ ■ • denote oscillator quantum numbers, 
s z - +i (s z = — i) the projection of electron spin f (J.) Ez = 
g*HBBs z the spin Zeeman energy, and u> ± - a)/,+w c /2 is defined 
in terms of the hybridized frequency w/, = (oj^ + a> 2 /4) 1 ^ 2 and 
the cyclotron frequency a> c — eB/m*. . 

The eigenstate \n,m, q) possesses the orbital angular mo- 
mentum l z = h(n - m) and the parity P = 1 (P = -1) 
with respect to z-axis for the even (odd) q number. The 
wave function of the lowest orbital is given by ifoooO") = 
exp [- ((x 2 + y 2 ) /l 2 h + z 2 ll\) /4] / (2n 3/ \ *J%) with the charac- 
teristic lengths of the wave function extents If, = y/h/2m*LOi, 
and l z = sjh/2m*a> z , from which one can generate the wave 
functions of any other excited states by successively applying 
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Figure 1 : Single-electron energy spectrum as a function of aspect ratio a with 
no magnetic field. The dot diameter L v =50 nm and the dot height L z varies 
from 40nm to 150nm. The finite difference results are shown in the inset for 
comparison. 



Figure 2: Schematic diagram of the lowest two orbitals occupied by two rel- 
evant electrons for the cases: (a) isotropic (a = 1) and B = 0; (b) "rod- 
like" (a > 1) and B = 0; (c) with nonzero electronic g-factor and B ± 0. 
2Ez = g'/isB indicates the spin Zeeman energy splitting. 



raising operators Q5D - For most synthesized NWQDs, the di- 
ameter of the cross section is of the scales 50 nm, while the 
length of wire could be tunable over a wide range from 10 nm 
to 300 nm jf?]. To characterize the geometry of NWQDs, we 
define the aspect ratio parameter, 



(3) 



according to the extents of the wave function. 

Figure Q] shows the single electron energy spectra, as a func- 
tion of a, of NWQDs at zero magnetic field. To examine the 
validity of the model, we carry out a numerical finite differ- 
ence (FD) simulation for the electronic structure of InAs/InP 
heterostructure NWQD, as shown in the inset of Figure [T] 
The InAs NWQD is embedded in InP barriers with the diam- 
eter L x = 50 nm and varying the length from L z = 40 nm to 
L z = 150 nm. The effective mass m* = 0.023mo and the bar- 
rier offset Vi, = 0.6 eV are taken |7|. The confining strength 
parameter Hljq is fit by the ground state energy from the FD 
simulation. We set Ha>o = 13.3 meV for L x - 50 nm NWQDs 
throughout this paper. The numerically calculated energy spec- 
trum is in good agreement with that given by parabolic model. 
The schematic illustration of the engineered single electron en- 
ergy levels of NWQDs by tunable a and applied B is shown in 
Figure [2] 

2.2. Interacting NWQD 

To investigate the few-electron interaction effects in a 
NWQD, we express the few electron Hamiltonian in second 
quantization as 



W(r L in- 
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i,cr 
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(4) 



i jkl,crcr' 



where ;, j, k, I stand for the composite indices of single elec- 
tron orbitals (e.g. |/) = |«,-, m,-, qi)), cr =1 / I denotes the 
electron spin with s z = +1/ — 4, and c\ (c JO -) is the electron 



creation (annihilation) operators. The first (second) term on 
the right hand side of Eq.© represents the kinetic energy 
of electrons (the Coulomb interactions between electrons) 
and the Coulomb matrix elements are defined as (ij\V\kl) = 
e 1 (AitK)- 1 f J rfridr 2 ^(n)iA*(r 2 ) (|n - r^y 1 ^(r 2 )^(n), 
where k is the dielectric constant of dot material ( k = 15.15enis 
taken for InAs throughout this work). After lengthy derivation, 
for NWQDs with a > 1, we obtain the following formulation 
of the Coulomb matrix elements: 



(nitniqi; njmjqj\V\n0nkqk\ nflniqi) 
1 

(-r)- 
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where we have defined u = m,- + mj + m + n^ - (p i +p 2 +P4 +ps), 
v = (q t + qi + qj + q k )- 2(p 3 + p 6 ), R L = (m ; + m j) - (n t +nj) = 
-(L Zii + L zJ ), R R = (mi + m k ) - (tii + n k ) = -(L zJ + L Zik ), 
x = u) z la>h, and iF\ is me hypergeometric function. The 5- 
functions d qi+q . +qi+qk>even and (5« t ,R fi in the formulation ensure 
the conservation of the parity with respect to z-axis and the z- 
component of angular momentum of system L z , respectively. 
The formulation of Eq. (0 is reexamined by computing the 
Coulomb integral numerically. 
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Figure 3: Correlated two-electron energy spectrum as a function of magnetic 
field in a NWQD with diameter L x = 50 nm and aspect ratio a = 3. 

2.3. Exact diagonalization 

The energy spectrum of an interacting two-electron NWQD 
is calculated using the standard numerical exact diagonalization 
technique [8]. The numerically exact results are obtained by in- 
creasing the numbers of chosen single electron orbital basis and 
the corresponding two-electron configurations until a numerical 
convergence is achieved. In our full configuration interaction 
(FCI) calculations, we chose the typical orbital number from 20 
to 26 and the number of corresponding configurations is from 
190 to 325. In order to highlight the Coulomb correlations, we 
also carry out the partial CI (PCI) calculations in which only 
the lowest energy N e configuration is taken and compare the 
PCI results with those obtained from the FCI calculations. 



3. Results and discussion 

Figure |3]presents the FCI result of magneto-energy spectrum 
of two interacting electrons in a NWQD with a — 3. The ST 
transition of the two-electron ground state is shown to happen 
as Z?st ~ 0.9 T. As the applied magnetic field is weak, the spin 
Zeeman splitting is small and the two electrons mostly dou- 
bly fill the lowest S-orbital. With increasing magnetic field in- 
creases, the energy difference between triplet and singlet states 
of the two electrons decreases because of increasing spin Zee- 
man and exchange energies, both of which energetically fa- 
vor the triplet states \T + ). As the applied field is higher than 
B c ~ 0.9T, the ground state of two electrons transit from the 
singlet state to the triplet one. 

Now we turn to study the spin singlet-triplet splitting as a 
function of a, defined by A ST , = Ep - E$, with i = -, 0, 
and + corresponding to the T~, T°, and T + triplet states, re- 
spectively. In Figure EJa), we show the A ST , as a function of 
aspect ratio a under a fixed magnetic field B — 1 T. In the non- 
interacting case, Ag Ti are shown to decrease monotonically with 
increasing aspect ratio a. Since only T + energy is decreased by 
spin-Zeeman term, the ST transition could only occur between 
S and T + states. That is, only A sx + crosses zero as a is very 
large, while A ST o remains positive always. Thus below we shall 
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Figure 4: (a) Singlet-triplet splitting A ST , as a function of aspect ratio a with 
B = 1 T. (b) The spin phase diagram for electrons making singlet-triplet transi- 
tion with respect to magnetic field B versus aspect ratio a. 

only consider A sx + for the discussion of ST transition. The non- 
interacting ST splitting can be derived as Ag T+ = ha>o/a 2 - 2Ez, 
explicitly showing the quadratic decrease of Ag T+ with respect 
to a. Accordingly, in the non-interacting picture, the critical 
aspect ratio ast where the ST transition occurs is predicted as 

However, the PCI calculation predicts a much smaller value 
of critical aspect ratio «st = 2.5. In the PCI result, the ST 
splitting is substantially reduced by the energy reduction of the 
T state due to the reduced direct Coulomb interaction and the 
negative exchange interaction between the two electrons in the 
state. The FCI calculation shows «st = 2.9, as indicated by 
the dashed vertical line in Figure Ufa). In fact, the difference 
in the values of Ast+ obtained from the PCI and FCI calcula- 
tions increases as a increases. This indicates that the Coulomb 
correlation effect tends to increase the ST splitting again and 
becomes even more pronounced in long NWQD with high a. 

Figure HJb) shows the calculated spin phase diagram of two- 
electron NWQDs with respect to the aspect ratio a and applied 
magnetic field B. The spin singlet and triplet phases, appearing 
in the low a-B and high a-B regimes, respectively, are distin- 
guished by the curve of B$j which show a monotonic decrease 
with a. For noninteracting electrons, the critical magnetic field 
can be derived as Z?st = fi&>o/gMBa 2 , showing a quadratic decay 
with a. 

In comparison with the non-interacting cases, the PCI cal- 
culations obtain the Z?st that is significantly reduced and goes 
to zero for a > 2.9. In the one-configuration approxima- 
tion used in the PCI calculation, the ST splitting is given by 
A ST + ~ A° T+ + Afj - V| x , where A° T+ is the ST splitting in the 
non-interacting cases, A™ = V* 1 - VJ? ir < is the direct energy 
difference between the triplet and the singlet states, and V™ is 
the exchange energy between electrons in the T + state. Accord- 
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ingly, we obtain B ST = (hwo/a 2 + A^ - V^)/g*fx B - In the large 
aspect ratio regime, the negative and V™ reduce the E T + 
and Z?st = results for a > 3. However, the FCI calculation 
predict larger and always positive B ST . In fact, as increasing 
a, the relative strength of electronic Coulomb correlations in- 
creases because of reduced hco z and strong configuration inter- 
actions. Such a-engineered Coulomb correlations energetically 
favor the singlet two-electron states and result in the non-zero 
Bst in the high aspect ratio regime. 

4. Summary 

In conclusion, a configuration interaction (CI) theory is 
developed for studying the magneto-energy spectra and the 
singlet-triplet transitions of two-electron NWQDs with longi- 
tudinal magnetic field B and tunable aspect ratio a. For short 
NWQDs of low aspect ratio a < 3, the ST transition behav- 
iors are dominated by the spin Zeeman, Coulomb direct and 
exchange energies, and can be well studied by using PCI cal- 
culation. However, our FCI calculations show the increasing 
importance of Coulomb correlations in long NWQDs with in- 
creasing aspect ratio a over 3. The FCI calculation present the 
spin phase diagram of a two-electron NWQD which are highly 
dependent on a, and suggests the controllability of singlet or 
triplet spin states by changing the aspect ratio of NWQD. 
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